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Abstract 

A method is developed to search for air showers initiated by photons using data 
recorded by the surface detector of the Auger Observatory. The approach is based 
on observables sensitive to the longitudinal shower development, the signal risetime 
and the curvature of the shower front. Applying this method to the data, upper limits 
on the flux of photons of 3.8 x 10~ 3 , 2.5 x 10~ 3 , and 2.2 x 1CT 3 km~ 2 sr _1 yr _1 above 
10 19 eV, 2 x 10 19 eV, and 4 x 10 19 eV are derived, with corresponding limits on the 
fraction of photons being 2.0%, 5.1%, and 31% (all limits at 95% c.l.). These photon 
limits disfavor certain exotic models of sources of cosmic rays. The results also show 
that the approach adopted by the Auger Observatory to calibrate the shower energy 
is not strongly biased by a contamination from photons. 



1 Introduction 



The search for photons in the ultra-high energy (UHE) cosmic-ray flux has 
been stimulated by the observation of cosmic rays with energies exceeding 
Eqzk ~ 6xl0 19 eV [1,2,3,4,5,6]. If these particles are due to cosmologically 
distant sources, the flux spectrum is expected to steepen above this energy. In- 
triguingly, a flux spectrum with no apparent steepening above -Eqzk has been 
reported by the AG ASA Collaboration [7]. To account for this observation and 
to circumvent the theoretical challenge of explaining particle acceleration to 
such energies, models involving new physics have been proposed in which the 
cosmic rays are created at the observed energies at relatively close distances 
from the Earth. These "top-down" models [8,9] may involve super heavy dark 
matter (SHDM) [10,11,12], topological defects [13], or neutrino interactions 
with the relic neutrino background (Z-bursts) [14]. A common feature of these 
models is the prediction of a substantial photon flux at highest energies. 
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The Auger Collaboration has recently reported a measurement of the cosmic- 
ray spectrum from the Auger South site showing a flux suppression above 
Eqzk [15]. The Auger method is based on a large surface array to collect the 
required statistics and a fluorescence detector to calibrate the energy scale. Us- 
ing this "hybrid" approach, the energy reconstruction is largely independent of 
hadronic interaction parameters and, in case of nuclear primaries, of the pri- 
mary mass composition. However, as explained later, the energy assignment 
from surface arrays can be substantially altered in the case of primary pho- 
tons. This would affect the reconstructed primary spectrum if a non-negligible 
number of the highest-energy events, where data from the fluorescence tele- 
scopes are sparse due to their ~10% duty cylce, was actually due to photons 
(see also [16]). It is worthwhile to note that the acceptance of fluorescence 
detectors (as also applied in the HiRes experiment [5]) can be altered in the 
case of photon primaries [17,18,19]. 

UHE photons can also act as tracers of the GZK (Greisen-Zatsepin-Kuzmin) 
process [20] of resonant photopion production of nucleons off the cosmic mi- 
crowave background. The corresponding photon fluxes are sensitive to source 
features (type of primary, injection spectrum, distance to sources ...) and 
to propagation parameters (extragalactic radio backgrounds and magnetic 
fields) [9,21,22,23,24]. 

Thus, the search for primary photons remains an important subject for various 
reasons [25], particularly 

• to set significant limits to the possible contribution of top-down mechanisms 
to the primary cosmic-ray flux; 

• to search for GZK photons, to prove the GZK effect and constrain source 
and propagation models; 

• to establish the maximum photon fraction in the primary flux, for which 
the energy estimate in the surface array detector would be altered; 

• to obtain input to fundamental physics, for instance, to probe quantum 
gravity effects in the electromagnetic sector [26]. 

Showers initiated by UHE photons develop differently from showers induced by 
nuclear primaries. Particularly, observables related to the development stage 
or "age" of a shower (such as the depth of shower maximum A" max ) and to 
the content of shower muons provide good sensitivity to identify primary pho- 
tons. Photon showers are expected to develop deeper in the atmosphere (larger 
X max ). This is connected to the smaller multiplicity in electromagnetic inter- 
actions compared to hadronic ones, such that a larger number of interactions 
is required to degrade the energy to the critical energy where the cascading 
process stops. Additionally, the LPM effect [27] results in a suppression of 
the pair production and bremsstrahlung cross-sections. Photon showers also 
contain fewer secondary muons, since photoproduction and direct muon pair 
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production are expected to play only a sub-dominant role. 

Searches for photons were previously conducted based on surface arrays [28,29,30,31,32], 
and limits to the fraction of photons were reported (see [25] for a review). The 
derivation of limits to the photon fraction using surface array data alone is 
an experimental and conceptual challenge (see also Section 2.3). Firstly, for 
conclusions on the fraction, the energy scales for photon and nuclear primaries 
are needed. These energy scales may differ from each other for surface arrays, 
and the difference between the scales may depend in a non-trivial way on 
primary parameters such as the shower zenith angle. Secondly, the energy re- 
construction of nuclear primaries suffers from substantial uncertainties due to 
our limited knowledge of high-energy hadron dynamics. 

Both issues can be resolved using the fluorescence technique, which is near- 
calorimetric and largely independent of simulating hadron interactions. A cor- 
responding approach has been developed and applied recently to obtain a 
first bound on the fraction of photons from data taken at the Auger Observa- 
tory [19]. 

In this work, using the larger number of events recorded by the surface array, 
we derive for the first time a direct limit to the flux of photons by searching for 
photon candidates and relating their number to the well-known exposure of the 
surface array. This avoids the need of simulating events inititated by nuclear 
primaries; only the photon energy scale is needed which can be simulated 
with much higher confidence. Two observables of the surface detectors are 
chosen which have significantly different behavior for nuclear primaries when 
compared to photons: the risetime of the recorded shower signal and the radius 
of curvature of the shower front. 



We also derive a limit to the fraction of photons. While the challenge of using 
two energy scales remains for this part of the analysis, hadron simulations can 
still be avoided by using the hybrid calibration [15] to reconstruct the energies 
of the observed events. 



The plan of the paper is as follows. In Section 2, the observables used in 
the analysis and their relationship with the composition of cosmic rays are 
explained. In Section 3, the simulation of UHE photons is considered. The 
method developed to distinguish events which are photon candidates using 
observables of the surface detector is detailed in Section 4. In Section 5, the 
results are presented. The conclusions are given in Section 6. 



9 



2 Observables 



The analysis in this paper is based on data taken during 21,400 hours of 
operation of the surface detector recorded in the period 1 January 2004 to 31 
December 2006. The surface detector, when completed, will have 1600 water 
Cherenkov detectors spaced 1.5 km apart and covering ~3000 km 2 [33,34]. 
Each water Cherenkov detector, or station, is a cylinder 1.2 m in height and 
3.6 m in diameter. Each detector is lined with a reflective container that holds 
12 tonnes of purified water and is fitted with three nine-inch photomultiplier 
tubes (PMTs) looking down into the water. 

When a relativistic particle passes through a station, Cherenkov radiation 
is emitted. The radiated photons then propagate through the water, being 
reflected at the station walls, and are either eventually absorbed or detected 
by a PMT. The signals from the PMTs are digitised by a flash analog to 
digital converter (FADC) which samples the signal every 25 ns. These digitised 
signals are then transmitted to a central data acquisition system where event 
triggers are built. Each event, then, has a detailed time profile Sj(rj,i) of 
the energy deposited in each station i at distance in the shower plane. 
The function s(r, t) depends in a complex way both on the parameters of 
the primary particle (energy, type, direction) and on the detector response 
to different secondary particles (particularly the electromagnetic and muonic 
shower components). 

In this work, we extract two relatively simple but robust observables from these 
data, noting that the wealth of information contained in the time profiles can 
further be exploited in future work. The observables, the radius of curvature 
of the shower front and the risetime at 1000 m core distance, were found 
to provide good discrimination between photon and nuclear primaries (see 
e.g. also Ref. [35]). In addition to the quantitative studies of these observables 
by means of the simulation-reconstruction chain, we will also sketch (in a 
simplified way) why these observables are indeed expected to differ between 
nuclear and photon primaries. 



2. 1 Radius of Curvature 



Due to geometrical reasons, the arrival of the first particles at lateral distance 
r from the axis is expected to be delayed with respect to an (imaginary) planar 
shower front (see also Fig. 1, left plot). For a particle that is due to an earlier 
interaction at height H along the shower axis and observed at r, the delay 
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from the longer path length can be approximated as 



1 , r 'i 

t — -(s/H 2 + r 2 — H) oc— (r«if). 



(1) 



The delay increases (for r -C H about quadratically) with r. Importantly, the 
delay decreases with increasing height H. Air showers with the first ground 
particles coming from relatively large heights will have smaller delays t at fixed 
distance r compared to showers where the registered particles originated from 
smaller heights. Compared to primary photons, showers from nuclear primaries 
develop higher in the atmosphere (smaller -X" max ) . Additionally, shower muons 
(much more abundant in showers from nuclear primaries) can reach the ground 
from still higher altitudes further reducing the time delay. Thus, for nuclear 
primaries smaller delays are expected compared to photon primaries. 

We make use of this relation by fitting a shower front (abstract surface with 
convex curvature defined by the fastest shower particles) to the measured 
trigger times Ufa) of the first particles registered at distances rj. In the present 
study, the shape of the shower front is approximated using a spherical model 
(in accord with Eq. (1)), and the radius of curvature R of the shower front is 
obtained by minimizing x 2 m the function 



where U is the trigger time for station % as defined in [36], to is the time of the 
shower in the center of curvature, a is the unit vector along the shower axis, 
Xi is the location of the station on the ground relative to the shower core, and 
at is the uncertainty in the shower arrival time [37]. In the determination of 
U, a software filter is applied to reduce contributions from spurious signals not 
related to the actual shower. 

2.2 Risetime 

Also the spread in time of the signal Sj(rj,£) registered at distance r^, which 
corresponds to the thickness of the local shower disk, can be extracted. Using 
Eq. (1), the difference of arrival times of particles originating from a height 
interval [Hi, Hi — AH] follows as 



x 2 = E 



[cfa - t ) - \Ra-Xi\} 2 



(2) 



At(H u AH)ocr 2 



1 1 



) 



r 2 AH 



Hi - AH H 



i 



Hi(Hi-AH) 



< At(H 2 ,AH) for H 2 < H L . 



(3) 
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Fig. 1. Illustration of geometrical effects on radius of curvature and risetime of the 
shower front. (Left) With respect to an imaginary planar shower front, particles 
arrive more delayed at distance r when originating from a smaller height H2 < H\. 
Correspondingly, the radius of curvature of the actual shower front is smaller in 
case of the deep developing photon primaries. (Right) The spread of arrival times of 
particles produced over a pathlength AH and arriving at distance r increases for a 
smaller production height H2 < H\. Correspondingly, the risetime of the shower is 
increased in case of the deep developing photon primaries. 

The spread of arrival times of these particles at fixed core distance increases 
for smaller production heights (see also Fig. 1, right plot). Accordingly, a 
larger spread is expected in case of the deep developing photon primaries 
(larger X max ). We note that in general, the situation is more complex. The 
time spread may depend on details of the previous shower development, par- 
ticularly also on the competition between the signals from the electromagnetic 
and muonic shower components which will be commented on below. Still, geo- 
metrical effects are essential in the relation between time spread and primary 
composition. 

In this study, we use the risetime £i/ 2 (1000) of the shower signal reconstructed 
for 1000 m distance and located along the line given by the projection of the 
shower axis onto the ground. First, the risetime t™ c 2 as (rj) of a single station 
is defined as the time it takes to increase from 10% to 50% of the total sig- 
nal deposited in that station. According to Eq. (3), for non- vertical showers 
a (moderate) dependence of t™ c 2 as (rj) on the internal azimuth angle of the 
stations within the shower plane is expected. This is because the height H 
measured along the shower axis is larger for those stations on the exterior side 
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of the shower compared to those on the interior side of the shower. To account 
for this, the observed ty^fc) are corrected depending on the internal azimuth 
angle C of that station: 

^/ r 2 (n) = C/r(n)-^-cosC (4) 

g = -66.61 + 95.13 • sec 9 - 30.73 • sec 2 9 + [0.001993 • sec 9 
-0.001259 • sec 2 9 + 0.0002546 • sec 3 9 - 0.0009721] • rf 

where the parameter g depends on distance r and primary zenith angle 9 
and is parameterised from the data, and ( is the clockwise angle between the 
projection of the shower axis on the ground and the line connecting the shower 
impact point and the station. 

It is also expected from Eq. (3) that the values tffoiji) depend on the distance 
Ti of the stations. We obtain the final risetime ti/ 2 (1000) of the shower by 
performing a fit to t c °j 2 {r^) using the function 

t 1/2 (r) = (40 + ar + br 2 ) ns . (5) 

The parameters a and b are determined for each event by fitting the station 
data (typical values are 50 ns km -1 and 100 ns km -2 respectively). The func- 
tion is anchored at 40 ns at r=0 as that is the mean single particle response 
in the water Cherenkov detectors. 

While geometrical effects connected to the different shower developments from 
nuclear and photon primaries are a main reason for the risetime difference 
(larger ti/ 2 (1000) in photon showers), again this sensitivity to composition is 
further strengthened by shower mouns which are more abundant in the case 
of nuclear primaries and can dominate the registered signal at larger zenith 
angles. As muons tend to arrive within a shorter time window compared to 
the electromagnetic component which suffers from multiple scattering, this 
further reduces the risetime ti/ 2 (1000) for nuclear primaries. 

2. 3 Energy 

As an energy estimator, the time-integrated energy deposit 5(1000) at 1000 m 
core distance is used [38]. However, for the same initial energy and direction 
the average 5(1000) from primary photons can be a factor >2 below that 
from nuclear primaries [39,40]. Reasons are the (typically factor ~4) smaller 
number of muons and, due to the later development, the steeper ground lateral 
distribution in primary photon showers. For a limit to the fraction of primary 
photons, the energy scales (transformation from 5(1000) to primary energy) 
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for both photon and nuclear primaries are required, while the determination 
of a limit to the flux can rely on the photon energy scale alone. 



The energy scale for nuclear primaries is based on the fluorescence technique 
by using events that are detected with both the surface detector and the fluo- 
rescence telescopes [41]. The energy scale for photon primaries (which induce 
almost purely electromagnetic cascades) is taken from simulations. Thus, both 
approaches are largely independent from assumptions about hadron interac- 
tions at high energy. 

Using a direct relationship between S'(IOOO) and primary energy for the photon 
energy scale results in a (relatively poor) resolution of about 40%. To improve 
this, a unique energy conversion for photons is applied that is described in 
detail in Ref. [40]. It is based on the universality of shower development [42], 
i.e. the electromagnetic part of the shower is expected to develop in a well- 
predictable manner for depths exceeding X max . In brief, for given values of 
,S(1000) and X max , the primary energy is estimated by 

g<lM = 1 .4( 1 + AX - 100 ) [1 + ( AJf - 1 ' ) V]- ( 6) 
£ 7 V 1000 ;[ V 340 ; J K ' 

with AX = X ground — X max , 



where .5(1000) is measured in units of vertical equivalent muons (VEM) [36], 
the photon energy E 1 is in EeV, and AX is in g cm -2 . Since X max is not directly 
measured by the surface detector alone, an iterative approach using Eq. (6) is 
taken to estimate the energy. After an initial guess of the photon energy using 
5(1000) alone, the typical X max of the photon showers at this energy is taken 
from simulations. With this estimate of X max , a new estimate of the photon 
energy is obtained using Eq. (6), and the procedure is repeated. The energy 
estimate is found stable after few iterations and an energy resolution of ~25% 
is achieved [40]. We use this improved estimation of the photon energy, but 
note that the main conclusions remain valid also when using a direct energy 
estimation. 



3 Monte Carlo Simulations 



The QED processes of LPM effect [27] and geomagnetic cascading ([43,35] and 
references therein) need to be considered for photon showers at highest energy. 
As mentioned before, the LPM effect leads to a suppression of the pair produc- 
tion and bremsstrahlung cross-sections and, thus, additionally increases the 
separation of photon and nuclear primaries in terms of X max (for a review of the 
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LPM effect and experimental observations of the LPM suppression, see [44]) [jj 
In case of geomagnetic cascading of UHE photons, the initial conversion of the 
UHE photon into an electron-positron pair can induce a "preshower" (mostly 
synchrotron photons plus electron-positron pair(s)) outside the atmosphere. 
The subsequent air showers from such "converted" photons develop higher in 
the atmosphere (smaller -X" max ) than air showers directly initiated by UHE 
photons do. As geomagnetic cascading becomes important at energies above 
~50 EeV at the southern site of the Auger Observatory, this process is of 
minor relevance for the bulk of data used in this analysis. 

The shower simulations were generated with the Aires simulation package 
(v2.8), which includes the LPM effect and geomagnetic cascading [45]. QGS- 
JET 01 [46] was used as the hadronic interaction model. The simulation of 
the water Cherenkov detectors uses the GEANT4 [47] simulation package 
along with specific code that handles PMT response and data acquisition 
electronics. The result is that the output of a simulated event is in a format 
that is identical to the data format recorded with the Auger Observatory. The 
shower reconstruction procedure used is the same for real events as it is for 
simulated events to avoid systematic differences at the reconstruction stage. 



4 Method 



In brief, the limit to the photon flux is obtained as follows. Selection cuts are 
applied to the data (and simulations) to ensure events of good reconstruction 
quality and a high acceptance of the detector to photons. Based on S'(IOOO), 
showers above a minimum primary energy are selected. This data set is then 
searched for photon candidates using ti/ 2 (1000) and R (see Section 2 for defini- 
tions). Simulations assuming photons are used to determine the corresponding 
selection efficiencies. From the number of photon candidates, the efficiencies 
with respect to photons, and the experimental exposure (obtained from the 
geometrical acceptance known from detector monitoring), the upper limit to 
the photon flux is derived. 

The criteria to select events of good quality are: 

• the station with the largest signal is surrounded by 6 active stations; 

• >5 stations used in the fitting of the lateral distribution function [48] out 
of which >4 stations have a non-saturated signal of >10 VEM (vertical 
equivalent muons) [36]; 



Even when artificially switching off the LPM effect, photon showers still have 
a significantly larger X max than nuclear primaries (differences >150 g cm 2 above 
10 19 eV) and a smaller number of muons. 
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Fig. 2. Photon detection and reconstruction efficiency (right hand scale) as a function 
of the energy (in EeV) and zenith angle of the primary photon. The analysis is 
restricted to a minimum energy of 10 EeV and zenith angles greater than 30° and 
less than 60° (0.866 > cosfl > 0.5). 

• reduced x 2 < 10 (x 2 from Eq. (2)). 

The first cut restricts the analysis to well-contained events, eliminating in 
particular events near the border of the array. It affects the geometrical ac- 
ceptance only. The multiplicity criterion in the second cut is important also 
to ensure a good reconstruction of ti/2(1000) and R. As the multiplicity is 
related to primary energy, this cut also affects the energy-dependent accep- 
tance of the array to photons. The third cut rejects the extreme tail of the x 2 
distribution when reconstructing R, removing ~4% of data. As noted before, 
the assumption of a spherical model used in Eq. (2) is a simplification and, 
thus, not expected to provide a perfect description of the complex features of 
the shower front. This cut restricts the analysis to events where a single value 
of R can be reasonably extracted. It has been checked with simulations that 
no bias to photons is introduced this way. 

As can be seen from Fig. 2, the resulting photon efficiency drops to small 
values below ~10 EeV. At higher energy, near-vertical photons can also fail 
the station multiplicity cut due to their deep development. Therefore, the 
analysis is restricted to 

• primary energies >10 EeV; 

• primary zenith angles of 30—60°. 
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Fig. 3. Parameterization of the mean behavior of R and tx/2 for 20 EeV primary 
photons as a function of the zenith angle using QGSJET 01 [46] or SIBYLL 2.1 [49]. 
The rms values are indicated for the case of QGSJET 01. An increase (a decrease) 
of R (of ti/2) with zenith angle is qualitatively expected from Eqs. (1) and (3) due 
to the generally longer path lengths to ground in case of larger inclination. Real 
events of 19-21 EeV (photon energy scale) are added. The significant deviation of 
the observed values from those expected for primary photons is visible. 

Events with zenith angles below 60° are selected here since inclined show- 
ers require dedicated algorithms for an optimum reconstruction [50] (this cut 
might be relaxed in the future). 

The search for photon candidates makes use of £i/ 2 (1000) and R and consists 
of the following steps. Firstly, the deviation A x of the observable x (with 
x = ti/2 or R referring to risetime or radius of curvature, respectively) from 
the mean value x 7 predicted for photons is derived in units of the spread <j x ,j 
of the observable x, 

g-g 7 (g(1000),g) 
^(5(1000), 0) ' 1 ' 



where ^(S^IOOO), 9) and cr a;]7 (<S'(1000), 9) are parameterized from simulations 
using primary photons. In Fig. 3, examples are shown for these parameteriza- 
tions of the observables along with distributions of real events. 

Secondly, we combine the information contained in the quantities A tl/2 and 
Ar by performing a principal component analysis [51], leaving a more sophisti- 
cated statistical analysis for the future. To determine the principal component 
(defined as the axis with the largest variance), 5% of the real events are used 
together with results from photon simulations, see Fig. 4. For the simulations, 
a power law spectrum of index -2.0 has been assumed (see below for other 
indices). The remaining 95% of the data are then projected onto the principal 
axis along with the simulated photons. 

This procedure allows the a priori definition of a simple cut in the projected 
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Fig. 4. The deviation from a photon prediction for 5% of the data (closed squares) 
and simulated photon events (crosses). The solid grey line is the principal component 
axis identified using the limited set of real showers while the dashed line is the 
axis perpendicular to the principal component. The minimum energy is 10 EeV 
(£ 7 > 10 EeV). 



distribution to finally obtain photon candidate events. The cut was chosen 
at the mean of the distribution for photons, such that the efficiency of this 
cut is f — 0.5 by construction. Any real event falling above this cut will 
be considered a photon candidate. We note that such photon candidates, if 
occuring, can not yet be considered as being photons, as they actually might be 
due to background events from nuclear primaries. A presence of background 
events would result in weaker upper limits (larger numerical values) in the 
analysis approach adopted here. 

Finally, an upper limit on the number of photons Af^ h at confidence level CL is 
calculated from the number of photon candidate events iV 7 above a minimum 
energy, E min . The upper limit on the flux or fraction of photons above a given 
energy is based on Af^ h along with the integrated efficiency e of accepting 
photons, the photon selection cut efficiency (/ = 0.5), and either the exposure 
A of the detector for the flux limit: 

* CL (S >£„,„)= \ XA > ' (8) 
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or the number of non-photon candidate events iV non _ 7 in the data set for the 
fraction limit: 

A/; CL (£ 7 > E mm ) x I x \ 
Tcl{E> E min ) = — — - — — — — — r . (9) 

7\ 7 J - J mm) 1~ J 'non— 7 V-^non— 7 -^minj 



In Eq. (8), the factor 0.95 is from the fact that only 95% of the data are used 
to determine the number of photon candidate events. The energy is labeled 
as either the energy according to the photon energy reconstruction, E 1 , or 
(required in Eq. (9)) the energy according to the non-photon energy recon- 
struction, i? non _ 7 . 

Experimentally, the limit $cl to the flux is more robust than the limit JF CL 
to the fraction due to the different denominators of Eqs. (8) and (9). For J-'cl, 
two energy scales are required; also, with increasing energy, the statistical 
uncertainty of the quantity (iV 7 + iV non _ 7 ) becomes large. For $cl, in contrast, 
the aperture is known to good (~3%) accuracy. 

Though the present work does not aim at extracting a composition of nuclear 
primaries, it is interesting to check whether the principal component axis found 
from real data and the separation along it reflects what would be expected 
if the bulk of the real data is due to nuclear primaries. In Fig. 5, the same 
simulated photon events are used as in Fig. 4 but the 5% of real data are 
replaced with a set of ~750 Monte Carlo proton and iron showers with an 
energy of 10 EeV. The separation observed in real data is both in the same 
direction and of a similar magnitude as that expected from simulated nuclear 
primaries. 



5 Results 



The data from 2004-2006 are analysed as described in the preceding section. 
The integrated aperture of the Observatory is 3130 km 2 sr yr for the angular 
coverage regarded in this analysis. Above 10, 20, and 40 EeV, for the energy 
scale of photons (in brackets for nuclear primaries), the data set consists of 
2761 (570), 1329 (145), and 372 (21) events. The measured values of i 1/2 (1000) 
and R are used to determine the projection on the principal axis. A scatter 
plot of this quantity vs. the primary energy is shown in Fig. 6, while in Fig. 7 
the corresponding distributions are plotted for the three threshold energies. 
No event passes the photon candidate cut. The upper limits on the photon 
flux above 10, 20, and 40 EeV are then 3.8 x 10" 3 , 2.5 x 10" 3 , and 2.2 x 
10~ 3 kmT 2 sr -1 yr -1 (at 95% CL). The limits on the photon fraction are 2.0%, 
5.1%, and 31% (at 95% CL) above 10, 20, and 40 EeV. In Tab. 1, all relevant 
quantities (number of events, efficiencies, resulting limits) are summarized. 
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Fig. 5. The black crosses are simulated photon showers while the squares are a 
mixture of Monte Carlo proton and iron with an energy of 10 EeV. For comparison, 
the lines shown in Fig. 4 (principal component axes) are added. The distribution of 
simulated nuclear primaries is similar to the distribution of real data seen in Fig. 4. 
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Table 1 

Results of the analysis searching for photon candidate events. The fraction and 
flux limits are integral limits above E m \ n (EeV), e is the efficiency of detection and 
reconstruction, $o.95 is in units of km -2 sr _1 yr , and all results are 95% confidence 
level. 

From Fig. 6 it can also be seen that the separation of data and photon pri- 
maries increases with energy. In particular at highest energies above Eqzk 
for the photon energy scale, there is no indication for photon-initiated events. 
Thus, the absence of photons, within the improved limits placed in this work, 
shows that the method applied by the Auger Observatory to calibrate the 
shower energy is not strongly biased by a photon "contamination". 



We studied potential sources of systematic effects in the analysis. To determine 
the efficiency to photons and to establish the photon candidate cut, a primary 
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Fig. 6. The deviation of data (black crosses) and photons (open red circles) from 
the principal component as a function of the primary energy (photon energy scale). 
Data lying above the dashed line, which indicates the mean of the distribution for 
photons, are taken as photon candidates. No event meets this requirement. 

photon spectrum of power law index -2.0 has been used in the simulations, 
motivated by predictions from top-down models in (e.g. in Ref. [10]). The effect 
of changing the power law index to -1.7, -2.5, and -3.0 has been investigated. 
The number of events which are photon candidates is unchanged (along with 
the number of non-photon candidate events), but the correction for the photon 
efficiency changes. Specifically, for a steeper input spectrum (increased fraction 
of lower-energy photons), the efficiency decreases. The summary of the results 
can be seen in Table 2. For 10 EeV threshold energy, limits change from 
(3.8->5.5) xlO -3 km~ 2 sr" 1 yr" 1 for the flux and from (2.0->2.9)% for the 
fraction. The differences get smaller with increased threshold energy. 

The photonuclear cross-section used in the simulation is based on the Particle 
Data Group (PDG) extrapolation [52]. For an increased cross-section, more 
energy would be transfered to the hadron (and muon) component which could 
diminish the separation power between data and primary photons [53]. From 
unitarity constraints, the cross-section is not expected to exceed the PDG ex- 
trapolation by more than ~75% at 10 EeV [54]; at 10 15 eV, where the difference 
in cross-section would have a greater impact on the shower development, the 
maximum difference is ~20%. From simulations with modified cross-sections 
it was verified that this leads to a negligible variation of the average values of 
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Fig. 7. Distribution of real events (closed squares) along with simulated photon 
events (open circles) for the projection on the principal component axis. The photon 
candidate cut is set at the mean of the distribution for photons and is shown as the 
dotted line. The plots are made requiring a minimum energy (according to the photon 
energy converter) of 10 EeV (top-left), 20 EeV (top-right), and 40 EeV (bottom). 
Distributions are normalised to unity at maximum. 

the discriminating variables used in the current analysis. 

The simulations have been performed with AIRES using the QGSJET had- 
ronic interaction model. As a cross-check, calculations with CORSIKA [55] / 
QGSJET and AIRES / SIBYLL were conducted, both of which show reason- 
able agreement to the AIRES / QGSJET case. As the cascade initiated by 
primary photons has an almost pure electromagnetic nature, indeed no sig- 
nificant effect is expected when changing to another interaction model. This 
minor dependence of the results on the details of hadron interactions, which 
are largely uncertain at high energy, is an important advantage of searches for 
primary photons. 

The new limits are compared to previous results and to theoretical predictions 
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Table 2 

Results when changing the exponent (a) in the power law of the simulated spectrum. 
The default value is 2.0. The efficiency of detection and reconstruction is on the left, 
the resulting limit on the fraction of photons is on the right, and the limit on the 
integrated flux is listed in the middle in units of km -2 sr _1 yr _1 (95% CL). 

in Fig. 8 for the photon flux and in Fig. 9 for the photon fraction. We placed the 
first direct limit to the flux of UHE photons (an earlier bound from AGASA, 
about an order of magnitude weaker than the current bounds, was derived 
indirectly via a limit to the fraction and the flux spectrum [29]). In terms of 
the photon fraction, the current bound at 10 EeV approaches the 10~ 2 level 
while previous bounds were at the 10 _1 level. 

A discovery of a substantial photon flux could have been interpreted as a signa- 
ture of top-down models. In turn, the experimental limits now put strong con- 
straints on these models. For instance, certain SHDM or TD models discussed 
in the literature (SHDM and TD from Ref. [211 based on the fragmentation 
calculations of Ref. [11], SHDM' from Ref. [12]_£j) predict fluxes that exceed 
the limits by a factor ~10. It should be noted that a simple rescaling of the 
flux predictions from top-down models, which were motivated by and based on 
the energy spectrum observed by AGASA, would reduce the predicted photon 
flux by only a factor ~2 which would still overshoot our experimental limit 
by a factor ~5 at 10 19 eV. While a minor contribution from top-down models 
to the observed UHE cosmic-ray flux might still be allowed within the limits 
derived in this work, current top-down models do not appear to provide an 
adequate explanation of the origin of the highest-energy cosmic rays (see also 
Ref. [56] for a comparison of photon flux predictions to the Auger limits for 
different top-down model parameters). 

In acceleration models, photon fluxes are usually expected to be a factor 2 or 
more below the current bounds (cf. the GZK photon predictions in the Figs. 8 
and 9 from Ref. [21]). Such fluxes can be tested with future data taken at the 
Auger Observatory (see also Ref. [25]). After five years of operation with the 
complete surface detector, sensitivities at the level of ~4x 10" 4 km -2 sr _1 yr -1 
for the integrated flux and ~0.7% for the fraction of photons above 20 EeV 

2 Two others of the eight photon flux spectra calculated in Ref. [12] from crypton 
decays may still be compatible with our limits within a factor ~2. 
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Fig. 8. The upper limits on the integral flux of photons derived in this work (black 
arrows) along with predictions from top-down models (SHDM, TD and ZB from 
Ref. [21], SHDM' from Ref. [12]) and with predictions of the GZK photon flux [21]. 
A flux limit derived indirectly by AG AS A ("A") is shown for comparison [29]. 

(95% CL) could be reached. 



6 Conclusions 

Using data from the surface detector we obtained 95% c.l. upper limits on the 
photon flux of 3.8 x 10~ 3 , 2.5 x 10~ 3 , and 2.2 x 10~ 3 km~ 2 sr _1 yr -1 above 
10 19 eV, 2 x 10 19 eV, and 4 x 10 19 eV. These are the first direct bounds on the 
flux of UHE photons. For the photon fraction, limits of 2.0%, 5.1%, and 31% 
were placed. 

These limit improve significantly upon bounds from previous experiments and 
put strong constraints on certain models of the origin of cosmic rays. Current 
top-down models such as the super-heavy dark matter scenario do not appear 
to provide an adequate explanation of the UHE cosmic rays. In bottom-up 
models of acceleration of nuclear primaries in astrophysical sources, the ex- 
pected photon fluxes are typically well below the current bounds. An astro- 
physical origin of UHE cosmic rays is also suggested by the recent discovery 
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Fig. 9. The upper limits on the fraction of photons in the integral cosmic-ray flux de- 
rived in this work (black arrows) along with previous experimental limits (HP: Hav- 
erah Park [28]; Al, A2: AGASA [29,30]; AY: AGASA- Yakutsk [31]; Y: Yakutsk [32]; 
FD: Auger hybrid limit [19]). Also shown are predictions from top-down models 
(SHDM, TD and ZB from Ref. [21], SHDM' from Ref. [12]) and predictions of the 
GZK photon fraction [21]. 

of a correlation of UHE cosmic rays with the directions of nearby AGNs [57]. 
Concerning the method of energy calibration as applied by the Auger Obser- 
vatory, the photon bounds derived in this work show that there is no strong 
bias due to a contamination from UHE photons. 

With the data accumulating over the next years, and particularly when com- 
plementing the Auger southern site by an extended northern one, the flux 
levels expected for GZK photons may be in reach. 
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